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Charmed hadron physics in quenched anisotropic lattice QCD* 
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We investigate the anisotropic lattice with O(a) improved quark action as a candidate of framework in which 
we can treat both the heavy and light quark region in the same manner and systematically reduce the systematic 
uncertainties. To examine applicability of anisotropic lattice, we calculate the charmed meson spectrum and decay 
constants in quenched approximation. We find consistent result with most advanced results on isotropic lattices. 



1. Introduction 

Recent experimental progress in flavor physics 
to look for the effect of new physics strongly 
requires precise theoretical predictions from the 
standard model. However, model independent 
calculation of hadronic matrix elements is difficult 
because of nonperturbative nature of QCD. The 
lattice QCD simulation is one of most promising 
approaches in which the systematic uncertainties 
can be reduced systematically jj]]. The ultimate 
goal of this work is to construct a framework for 
lattice calculations of hadronic matrix elements 
in a few percent systematic accuracy, as required 
by the experiments in progress. 

In a lattice calculation of heavy quarks such 
as the charm and bottom, one need to either 
avoid or control the large discretization error of 
0((mQa) n ). Although the effective theoretical 
approaches |||| successfully describe the matrix 
elements with heavy quarks within 10% accuracy, 
further improvement of accuracy is difficult due to 
perturbative error. The nonperturbative renor- 
malization method are not directly applicable 
to a determination of the coefficient of 0(a) im- 
proving term in the region of tuq ~ a -1 because 
of a strong mass dependence. Therefore a new 
framework is desired for which a systematic non- 
perturbative improvement can be performed. 

The anisotropic lattice, on which the temporal 
lattice spacing a T is finer than the spatial one a a , 
is a candidate of such a framework. A high tem- 



poral lattice cutoff may enable a relativistic treat- 
ment of heavy quark while keeping computational 
requirement tractable. The 0(a) improvement 
coefficients determined in the light quark region 
in good precision may directly be available for a 
heavy quark with itiq <C a^ 1 . Whether these 
observations practically hold should be examined 
numerically, as well as in the perturbation theory. 

In this study, we examine the applicability of 
anisotropic lattice approach numerically by com- 
puting the charmed meson masses and decay con- 
stants. Although presently achieved accuracy is 
far from our goal, the result is encouraging for 
further development in this direction. This pa- 
per shows our preliminary result without techni- 
cal details and evaluation of the systematic uncer- 
tainties, which are fully discussed in other publi- 
cations 1^-0]. 

2. Anisotropic lattice quark action 

The quark action we adopt has the following 
structure in the hopping parameter form J^,||]: 



S F = 5^(z) {*«,»- *r [ {l--n)Ui{x)8 
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where k ct and the spatial and temporal 

hopping parameters, and are related to the bare 
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quark mass and bare anisotropy as given below. 
r is the spatial Wilson parameter, and ce and cb 
are the coefficients of clover terms which remove 
0(a) errors. Although the explicit Lorentz sym- 
metry is removed due to the anisotropy, in prin- 
ciple it can be restored for physical observables 
up to errors of 0(a 2 ) by proper tuning of k ct /k t , 
r, ce and cb for a given K a . The action is con- 
structed in accord with the Fermilab approach 
H and hence applicable to an arbitrary quark 
mass. However, a mass dependent tuning of pa- 
rameters is difficult beyond perturbation theory. 
This may be circumvented by taking a~ x 3> tuq, 
with which the mass dependence of parameters 
are expected to be small, and it suffices with those 
values in the chiral limit. 

In present study, we vary only two parame- 
ters tic and k t with fixed other parameters. We 
set the Wilson parameter as r = l/£ and the 
clover coefficients as the tadpole-improved tree- 
level values, ce = ^-/u a u 2 , and cb = The 
tadpole improvement || is achieved by rescal- 
ing the link variable as Ui{x) — > Ui(x)/u a and 
U^{x) — > Ui(x)/u T , with the mean-field values of 
the spatial and temporal link variables, u„ and 
u T , respectively. Instead of k ct and k t , we intro- 
duce k and 7f as 

- = — 2( 7F + 3r-4) = 2(m 7f +4), 

7_F = K T U T jK a U a . (2) 

The former plays the same role as on the isotropic 
lattice, and the latter corresponds to the bare 
anisotropy. 

On an anisotropic lattice, one must tune the 
parameters so that the anisotropy of quark field, 
£f, equals that of the gauge field, £g : 

?f(/3,7g;k,7f) = €g(P, 1G\ k, jf) = £. (3) 

Although £q and £p are in general functions of 
both of gauge parameters (0, 7 g) and quark pa- 
rameters (k, 7 _f), on a quenched lattice one can 
determine £g — £ independently of k and 7 _f, and 
then tune jp so that a certain observable satis- 
fies the condition (||). In this work, we define 
£f through the relativistic dispersion relation of 
meson, 

E 2 (p) = m 2 +p 2 /eF+o[( P 2 n (4) 



for calibration. In the above expression, the en- 
ergy and mass E and m are in temporal lattice 
units while the momentum p is in spatial lattice 
units. £p converts the momentum into that in 
temporal lattice units. 

3. Numerical simulation 

Numerical simulations are performed on two 
quenched lattices of sizes 16 3 x 128 and 20 3 x 
160 with the Wilson plaquette action at = 
5.95 and 6.10, respectively, with the renormal- 
ized anisotropy ^ = 40. The values of bare 
anisotropy are chosen according to a numerical 
result performed in one percent accuracy in |lC(] . 
The lattice scales are set by K* meson mass, and 
result in a" 1 = 1.525(27) GeV and 1.817(22) GeV 
for (3 = 5.95 and 6.10, respectively. These val- 
ues deviate to about 10% from scales determined 
with other physical quantity, say hadronic radius 
ro, and represent the effect in neglecting dynam- 
ical quarks. The mean field values u T and u a are 
obtained in the Landau gauge || . 

The calibration was done in the quark mass re- 
gion below around charm quark mass in 0. It 
was found that the quark mass dependence of 
bare anisotropy 7^,, with which £p = £ holds, 
is actually small, and 7^ is well fitted to a linear 
form in m, 2 , 

^K) = Co + C 2 m^, m? = !(!-£). (5) 

The fit results in the values (Co,C2, K c) = 
(0.2490(8), 0.189(15), 0.12592(8)) at = 5.95, 
and (0.2479(9), 0.143(14), 0.12558(9)) at = 
6.10. These values are determined in 1% statisti- 
cal accuracy, while in the chiral limit additional 
1% error exists due to the form of extrapolation. 
It was also shown that the systematic discretiza- 
tion errors decrease toward the continuum limit. 
The light hadron spectra with the obtained pa- 
rameters are consistent with previous works on 
isotropic lattices. 

For heavy quark, we use four values of k cov- 
ering the physical charm quark mass with 7^ ac- 
cording to Eq. (||). As the light quark, we use 
three values of k with 7^ at the massless limit, 
which correspond to the lightest three values used 
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Table 1 

The heavy-light meson masses and decay con- 
stants for physical quark masses in physical units. 



[GcV] 


(3 = 5.95 


(3 = 6.10 


m D s 

m D » - m D 
m D . - m Ds 


1.9740(28) 
0.1018(85) 
0.0980(42) 


1.9760(27) 
0.1022(74) 
0.0921(39) 


U 

fK 
fo 

Id. 


0.1655(35) 
0.1868(29) 
0.2509(64) 
0.2863(40) 


0.1462(35) 
0.1665(30) 
0.2245(58) 
0.2570(37) 


f D/fa 

/d s / /d 


1.516(44) 
1.141(15) 


1.536(49) 
1.145(16) 



0.45 
0.40 
0.35 
0.30 
0.25 
0.20 



(3=6.10 (K* meson mass scale) 

o m,=0 
□ m.=m. 
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Figure 1. The heavy- light decay constant multi- 
plied by \pm in physical units. 



in H for light hadron spectroscopy. The meson 
correlator is calculated with the local meson op- 
erator. For the decay constant, we focus on the 
pseudoscalar mesons. While one needs the match- 
ing constants to convert the lattice result into the 
continuum theory, we only perform the tadpole- 
improved tree-level matching in this paper. 

We first extrapolate the result to the chiral 
limit, linearly in light pseudoscalar meson mass 
squared. The physical (u, d) and s quark masses 
are defined with the massless limit and the K me- 
son mass. Then they are fitted to the quadratic 
form in 1/toh, ttlh the PS heavy-light meson 
mass, and interpolated to the physical mjj me- 



son mass. The decay constant of heavy-light me- 
son is extrapolated and interpolated in the form 
of fa yJm-H ■ The result after chiral extrapola- 
tion is shown in Figure The gross feature of 
ran dependence is consistent with previous works 
q|. The masses and decay constants for physical 
quark masses are listed in Table 0. 

The hyperfine splittings are consistent with 
previous works, while 35% less than the exper- 
imental values. On the other hand, the splitting 
m D B — m D is close to the experimental value. The 
j3 dependence of decay constants are rather large, 
and considered as of the discretization error and 
renormalization effect, as well as the uncertainty 
in setting the scale on quenched lattice. Tak- 
ing ratios of decay constants, (3 dependence is 
largely cancelled. The result is consistent with 
the previous lattice works (!]]. For more quanti- 
tative discussion, we need to discuss more care- 
fully on the systematic uncertainties due to the 
anisotropy, and to include the matching of cur- 
rent with continuum theory. As our conclusion, 
the result of numerical simulation is encouraging 
to pursuit more quantitative analysis in this di- 
rection. 
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